Our goal is to determine light curves and period-luminosity 
(PL) relations in the near-IR, to assess the similarity between 
cluster and field pulsators, and to examine the predictive ca- 
pability of current pulsation models. The light curves are ob- 
tained from multiwavelength broadband J,H,K S photometry of 
Cepheids in both clusters, with periods previously established 
from optical photometry. Mean magnitudes for the Cepheids 
are used to construct PL relations in the near-IR. The proper- 
ties in the PL planes are compared with the behavior of field 
Cepheids in the LMC and with the predictions of recent pul- 
sational models, both canonical and overluminous. Cluster and 
field Cepheids are homogeneous and the inclusion of the clus- 
ter Cepheids in the field sample extends nicely the PL relation. 
The slope of the PL relation is constant over the whole period 
f^*) range and does not show - at least in the adopted IR bands 
< ^ - the break in slope at P ~ 10 d reported by some authors. A 
comparison with the predictions of pulsation models allows an 

■ estimate for the distance moduli of NGC 1866 and NGC 2031. 
/■""n \ The two clusters are found to lie at essentially the same dis- 

• tance. Fitting of theoretical models to the data gives, for the 
K filter, (m-M)o = 18.62+0.10 mag if canonical models are 
,— ' ■ used and (m-M)o = 18.42+0.10 mag if overluminous models 
| are used. On the basis of this result, some considerations on the 

■ relationship between the clusters and the internal structure of 
0^ the LMC are presented. 

: 
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I individual: NGC 1866 - globular clusters: individual: NGC 
\£) 2031 - Stars: distances 

o ' 

6 
& 

CO 



Astronomy & Astrophysics manuscript no. paper' final 


February 5, 2008 


(DOI: will be inserted by hand later) 





Infrared photometry of Cepheids in the LMC clusters NGC 1866 

and NGC 2031* 

Vincenzo Testa 1 , Marcella Marconi 2 , Ilaria Musella 2 , Vincenzo Ripepi 2 , Massimo Dall'Ora 2 , Francesco R. Ferraro 3 , 

Alessio Mucciarelli 3 , Mario Mateo 4 , and Patrick Cote 5 



INAF - Osservatorio Astronomico di Roma, Via Frascati 33, 00040 Monte Porzio Catone, ITALY 
INAF - Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, ITALY 
Dipartimento di Astronomia, Univ. di Bologna, Via Ranzani 1, 40127 Bologna, ITALY 
Department of Astronomy, University of Michigan, Ann Arbor, MI 48109-1090, USA 

Herzberg Institute of Astrophysics, National Research Council of Canada, 5071 West Saanich Road, Victoria, BC, V9E 2E7, 
Canada 



Received Juky 26, 2006/ Accepted October, 1 1 2006 

Abstract. Near infrared (IR) studies of Cepheid variables in the LMC take advantage of the reduced light curve amplitude 
and metallicity dependence at these wavelengths. This work presents such photometry for two young clusters known to contain 
sizeable Cepheid populations: NGC 1866 and NGC 2031. 



1. Introduction 

Cepheid variables are one of the cornerstones upon which the 
extragalactic distance scale is based. Since the pioneering work 
by Henrietta Leavitt in 1912 (Harvard College Observatory, 
Circ. No. 173), their importance has grown to the point that 
they are considered one of the most important distance in- 
dicators for galaxies within, and beyond, the Local Group. 
Recent HST studies have made use of Cepheids to measure 

the distances of galaxies beyond the Local G roup (e.g. the 

I II — 

Extragalactic Cepheid project Freedman et al. 2001) and to 

thereby calibrate a variety of secondary distance indicators. In 
this regard, the Magellanic Clouds (MCs) — and, in particular, 
the Large Magellanic Cloud (LMC) — offer a unique labora- 
tory for testing the Cepheid distance scale, being among the 
closest extragalactic environments with a resolved stellar pop- 
ulation. Since LMC Cepheids lie very nearly at the same dis- 
tance from the Sun, they offer an attractive way of sidestep- 
ping the problem of relative distance errors that plague stud- 
ies of the Period-Luminosity (PL) relation slope using Galactic 
samples. For these reasons, the LMC has long been an impor- 
tant first step in establishing the Cepheid distance ladder: i.e., 
the calibration of the period-luminosity (PL) relation in dis- 
tant galaxies often relies on an assumed LMC distance modu- 
lus (DM). For instance, in the Extragalactic Cepheid project, 
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a DM of (m - M) = 18.50 mag was adopted for the LMC 
iFreedman etaDfeOOll) . 

Among the most pressing issues to be considered when 
dealing with the Cepheid PL relation is its dependence on the 
metallicity, which has been questioned in a number of past 
works. From an observational perspect i ve, some authors (e.g . 
i Kennicutt et alJll998HSakai et alJEool iTammann etail fcoOS ; 
Ro maniello et all2005l) ar gue that such a dep endence is indeed 
present, while others (e.g. Gieren et al. 2005) suggest that any 
effect is small, perhaps i ll-defined (e.g. Gro enewegen & Sa laris 
2003: IStormetall2005l) . On the theoretical side, some authors 
again find a mild depende nce (if any) from linear nonadiabatic 
pulsa tion models (see e.g. lSaio & Gautsch\ll998HAlibert et alJ 
1999), whereas nonlinear convective pulsation models predict a 
significant depend ence on both metal content and helium abun- 
dance (see Fiorenti no et all2 002: Mar coni et al.l2005l and ref- 
erences therein). 

A determination of the PL relation presents different chal- 
lenges in different wavelength intervals. In the optical band- 
passes, the relation shows non-negligible scatter and the use 
of a color is required to reduce the dispersion (i.e., a Period- 
Luminosity-Color relation, or PLC). Moreover, the light curve 
amplitude is large so the measurement of the mean magnitude, 
on which both the PL and PLC relations are based, requires 
accurate photometry and well sampled light curves (although 
the large amplitude allows a better determination of the pe- 
riod because it is less sensitive to small scatter in the measure- 
ments). By contrast, the light curve amplitude is much smaller 
at near-IR wavelengths and the determination of the mean mag- 
nitude correspondingly easier. Moreover, the dependence on 
both the reddening and the color extension of the instability 
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strip is strongly reduced leading to a narrower PL relation than 
in the optical case. On the other hand, the reduced pulsation 
amplitude makes the period determination more difficult. In 
this work, in fact, we make use of periods already available 
from a variety of literature sources. 

Theoretical predictions based on nonlinear convective 
models confirm the great advantage of using near-IR filters to 
construct PL relations. As the wavelength increases from opti- 
cal to near-IR bands, PL relations become linear over the full 
period range spanned by the observed pulsators and their de- 
pendence on metallicity and the instrinsic width of the insta- 
bility strip decreas es significantly (see, e.g. lBono et alj fl999a: 
ICaputo et al.l200d) . 

In this context, classical Cepheids belonging to young stel- 
lar clusters play a particularly relevant role. In fact, not only 
they are at the same distance but are also characterized by the 
same age and chemical composition, thus offering a unique 
opportunity to investigate the uncertainties affecting both em- 
pirical approaches and theoretical scenarios. For the present 
study, we selected two young clusters having the largest known 
Cepheid population among LMC clusters: NGC 1866 and 
NGC2031. 

The first is the most massive young cluster in the age 
range ~ 100-200 Myr. It has been the subject of a very 
larg e number of papers in th e past, begi nning with those by 
lArp & Thackeravl i 19671) and iRobertsoiJ Jl974l) . Subsequent 
authors focused on studying the cluster either as a testbed of 
stellar evolution theory (e.g. iBrocato et alJ I 19891 IChiosi et alJ 
| l989HBrocato et alJl994HTesta et alJl999UWalker et alJ200ll 
Barmina et al. 2002; Brocato et al. 2003), as a Cepheid host 
JWelch et al|l99lUWelch & StetsoJl993tlGieren et alJll994t 
IWalkeil! 995UC,ieren et al J200ftlStorm et al l2004. or as a dy- 
namical laboratory IFischer et a J 1 19921) . It has also been the 
subject of a strong debate over the presence of convective over- 
shooting in intermediate-age stellar models, and on the frac- 
tion of binaries on the main sequence. Because the cluster lies 
in the outskirts of the LMC, field contamination is not severe. 
NGC 203 1 , on the other hand, is an order of magnitude smaller 
in mass and is located in the LMC bar, resulting in a high con- 
tamination of its stellar population. As an added complication, 
the surrounding field has an age comparable to the cluster (as 
is apparent in the CMD shown in Fig. Despite its large 
Cepheid population, it has not been extens ively studied in the 
past, except for an analysis of the CMD by Mo uld et"afl Jl993h 
who quote an age of ~ 140 Myr, very close to that of NGC 1 866. 

The metallicity of NGC 1866 has been estimated by various 
investi gators, with the spectroscopic determination of lHill etal] 
(2000K rFe/ Hl = -0.5 ± 0.1 con firming the photometric deter- 
mination of lHilker etaD dl995h . [Fe/H] = -0.46 + 0.18. The 
first — and, to our knowledge, only — spect roscopic de- 
termin ation of the metallicity of NGC 203 1 by iDirsch et alJ 
tOOd yielded [Fe/H] = -0.52 + 0.21, indistinguishable from 
NGC 1866. Moreover, both clusters have a fairly populated 
AGB that is clearly defined as a narrow sequence in the near-IR 
CMD (see Fig.Q, and the mo st luminous stars have been sub- 
ject of a spectroscopic study jMaceroni et alJ |2QQ2) that con- 
firms the evolutionary status of thermally-pulsing AGB stars at 
the tip of the AGB sequence. Taken together, these properties 



suggest that NGC 203 1 can be considered a low-mass analogue 
of NGC 1866. 

The paper is structured as follows: the observations and re- 
ductions are described in Sect. |2 Sect. [3] deals with the pre- 
sentations of the photometric sample of the Cepheids; in Sect. 
|4j IR light-curves and PL relations are presented along with a 
comparison to previous results and to the predictions of theo- 
retical models; we conclude with a summary of our results in 
Sect.El 

2. Observations, Reductions and Calibrations 

2.1. Observations 

Data were collected over four observing runs spanning an in- 
terval of nine years. Although the phase sampling is sparse and 
quite random, we show later that is not a serious issue in the 
near-IR. The log of observations is given in Tabled Despite 
the different telescopes and instrumentation used in this study, 
the observation strategy was the same for all the runs, following 
the standard procedure of taking dithered frames of the target to 
subtract the sky. The field size was 2'x2' for the first three runs, 
for which a 256 2 pixel NICMOS3 array was employed. When 
using the 1024 2 pixels camera SOFI at the ESO-NTT telescope 
(run #4), the field of view was 5'x5'. In Tabled columns record 
the telescopes and instruments used in each run, the total du- 
ration of the run, and the number of measurements for each 
cluster (last column: no. meas.). Note that, despite its five-night 
length, only a single point for the two clusters was obtained in 
run #4, as this was a su rvey program on LMC clusters (see. 
iMucciarelli et all J2006b ^. 

2.2. Reductions 

In order to check the relative consistency of all the data-sets, 
two widely used packages have been use d to generate o bject 
lists and magnitudes, namely D AOPHOT JStetsonl[l987l) and 
DoPHOT JSchechter et al.lll993h . The detailed reduction strat- 
egy was slightly different for the various runs. For runs #1 
and #2, sky flat-fields were taken during evening twilight; for 
run #3, dome-flats were taken using the ESO-MPI 2.2m tele- 
scope dome; and for run #4, dome-flats taken in the ESO- 
NTT dome were used together with screen-corrected images, 
as recommended in the SOFI manual. For all runs, dithered 
images were observed adopting off-source skies, taken some 
arcminutes away from the cluster center, and subtracting the 
sky frames to the single cluster images, before reconstructing 
the final scientific images by re-aligning the dithered source 
images. 

Bad pixel masks were used, whenever possible, by tak- 
ing advantage of the dithering technique to remove bad pixels 
from the final images. In practice, this technique did not al- 
ways work very well, especially for the LCO runs, when the 
chip had a fairly large number of dead pixels. However, these 
artifacts yield measurements that were easily identifiable from 
their large photometric errors. 

Photometry was performed with two different packages, 
in order to check for systematic differences due to the soft- 
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Run Id. Dates Site Tel. Instr. pix.sz field size no. meas. 

1 02-08 Dec 1992 LCO Swope lm IR camera 0.53" 136"x 136" 5 

2 12-16 Jan 1993 ESO MPI 2.2m IRAC2b 0.49" 125"x 125" 3 

3 01-08 Nov 1993 LCO Swope lm IR camera 0.53" 136"x 136" 5 

4 27-31 Dec 2001 ESO NTT SOFI 0.29" 295" x 295" 1 



Table 1 : Observing Log 



ware, a strategy adopted in the HST Cepheid distance scale 
projecl Freed man et alJ J200lh . It should be noted that these two 
packages differ fundamentally in approach: DAOPHOT carries 
out multi-profile PSF fitting on a group of objects whose PSFs 
interact with each other, while DoPHOT adopts a single-profile 
fitting starting from the brightest source of the list, then sub- 
tracting the measured object from the image and proceeding in 
order of decreasing brightness. 

When running DAOPHOT for all runs, all nights and all 
groups of three-filters (J, H, K^) observations, a master list 
of candidates was built by carefully aligning and summing 
the three images, then running daofind on the summed im- 
age. When running DoPHOT, the program was applied to the 
single images and the output lists cross-checked a-posteriori. 
However, within the internal fitting error, measurements were 
fully consistent for high-S/N sources. For fainter sources, dif- 
ferences were found to be somewhat larger, but always within 
3cr of the (larger) internal error. Calibration errors and aper- 
ture correction errors have been added to the internal errors to 
obtain the total photometric errors. 

2.3. Calibrations 

During runs #1, #2 and #3, standard stars were observed to ob- 
tain a primary calibration of the data. However, although the 
nights appeared to be photometric, a careful analysis of the 
calibrating standards revealed some inconsistencies that pre- 
vented us from obtaining a well calibrated sample. For this 
reason a further set of observations was taken during run #4, 
when conditions were photometric. In order to check the cali- 
bration, we cross-checked the output magnitudes against a 500 
stars from the 2MASS point- source catalog and with the pho- 
tometry of M ucciarelli et alJ d2006l) . based on the same data, 
and found the values to be consistent in both cases. The other 
runs were then transformed to that photometry of run #4 by 
applying a secondary calibration that, in all cases, consisted of 
a zero point term and a negligible color term. The magnitudes 
were then tran sformed to the LCO system, for consistency with 
Persson et alJ J2004I) . The relations used were obtained from 
CarpenterlfeOOll) : 



for each object, the list of measurements in all the images. In 
total, we obtained 14 photometric points, that were then av- 
eraged with a cr-clipping rejection algorithm to obtain the fi- 
nal photometry list of each cluster. Cluster Cepheids were then 
identified and their magnitudes extracted to construct the light 
curves. Image coordinates were transformed onto the system 
of run #4, for which an astrometric calibration was obtained 
from a subset of objects in the 2MASS catalog in common with 
those from our own catalog. The IR color-magnitude diagrams 
(CMDs) of NGC 186 6 and NGC 203 1 have be en presented and 
discussed elsewhere ( Mucciarelli et al.ll2006l). although a pre- 
liminary version was published in Maceroni et alJ ( 120021) . in or- 
der to identify the candidate TP-AGB stars for a spectroscopic 
study presented in that paper. Here, for the sake of complete- 
ness, we show the CMDs for the two clusters in Fig.^ in which 
our program Cepheids are identified by the filled squares. 
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Fig. 1: Near-IR color-magnitude diagrams for 
clusters. Program Cepheids are identified by the 
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the two LMC 
filled squares. 



Jlco = J2mass + 0.020 - 0.014 x (J 2M ass - K 2M ass) 

HlCO = H2MASS + 0.020 - 0.015 X (H2MASS _ K2mass) 
KlCO = K2MASS +0.015 

3. The Photometric Sample 

Photometry lists for all runs and all nights were cross-matched 
in order to obtain: (1) a master list of the photometry; and (2) 



3.1. The Cepheids 

The Cepheids for N GC 1866 were select e d from the lists 
of IWelch et ail dl99ll) and lWelch & Stetson! i 19931) . using the 
finding charts i n those papers . For NGC 2031, the sample was 
obtained from iMateol y 992) and from our own unpublished 
lists. Since our frames covered the central 2'x2' of the cluster 
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(with the exception of the SOFI field), the outermost Cepheids 
of NGC 1866 were not included in our IR list, but we iden- 
tified all the variables in the cluster core. Of the 15 Cepheids 
of our sample, 8 a re confirmed radial velocity members (see 
IWelch etal.lfl99ll) . The fact that the outermost variables of 
NGC 1866 are not part of the presented sample prevents com- 
p arisons with some o f the best-studied Cepheids, like the ones 
in St orm et all d2005l) — with which we have only a single ob- 
ject in common (HV 12202) — but the inclusion of the inner 
variables increases considerably the overall sample. We com- 
pared the J and K m agnitudes of HV 12202 with the values re- 
ported in Stor m et alJ (12005), finding a good agreement. Hence, 
rather than using single point photometry from run #4 for these 
objects, we preferred t o use the more precise mean values pub- 
lished bv lStorm et ail d2005t) . 

Because NGC 2031 is more compact, all the known ob- 
jects fall within our fields with the exception of variables V5 
and V7 that are not present in the current sample. Their peri- 
ods are, however, less precisely determined than the ones for 
NGC 1866. Indeed, as can be seen in Table[3] the periods have 
a typical precision of O.Old, which likely explains the overall 
poorer quality of NGC 2031 light curves. 



Name 


P (days) 


<J> 


o-J 


<H> 


o-H 


<K> 


o-K 


V4 


3.3157 


14.77 


0.05 


14.51 


0.03 


14.37 


0.04 


We8 


3.043 


14.86 


0.07 


14.67 


0.05 


14.51 


0.06 


V6 (ov) 


1.9433 


14.97 


0.04 


14.74 


0.05 


14.60 


0.05 


V8 (ov) 


2.0088 


14.93 


0.06 


14.68 


0.03 


14.55 


0.04 


V7 


3.453 


14.70 


0.06 


14.44 


0.04 


14.31 


0.04 


HV 12200 


2.7248 


14.96 


0.08 


14.71 


0.06 


14.53 


0.09 


HV 12202 


3.10112 


14.80 


0.09 


14.60 


0.06 


14.44 


0.04 


We2 


3.054 


14.87 


0.09 


14.58 


0.09 


14.41 


0.08 


We3 


3.045 


14.78 


0.06 


14.56 


0.08 


14.42 


0.10 


We4 


2.8604 


14.83 


0.13 


14.58 


0.13 


14.46 


0.12 


We5 


3.1763 


14.75 


0.08 


14.53 


0.05 


14.38 


0.10 


We6 


3.289 


14.76 


0.04 


14.52 


0.04 


14.37 


0.04 


WS5 


2.895 


14.87 


0.06 


14.61 


0.07 


14.50 


0.07 


WS9 


3.071 


14.64 


0.09 


14.49 


0.10 


14.30 


0.09 


WS11 


3.0544 


14.95 


0.10 


14.62 


0.17 


14.42 


0.16 



Table 2: Mean IR magnitudes for Cepheids in NGC 1866. 



Name 


P (days) 


<J> 


o-J 


<H> 


o-H 


<K> 


o-K 


VI 


3.07 


14.71 


0.08 


14.38 


0.07 


14.38 


0.06 


V2 


4.43 


14.43 


0.08 


14.03 


0.05 


13.98 


0.05 


V3 


3.96 


14.56 


0.09 


14.16 


0.08 


14.14 


0.09 


V4 


3.43 


14.73 


0.07 


14.37 


0.09 


14.39 


0.05 


V5 


3.32 


14.68 


0.12 


14.32 


0.12 


14.31 


0.12 


V6 


3.03 


14.85 


0.14 


14.50 


0.08 


14.51 


0.09 


V7 


3.13 


14.94 


0.08 


14.56 


0.08 


14.58 


0.09 


V8 


3.27 


14.75 


0.08 


14.39 


0.08 


14.36 


0.08 


V9 


2.95 


14.89 


0.07 


14.53 


0.08 


14.52 


0.06 


Vll 


2.82 


15.01 


0.12 


14.67 


0.07 


14.70 


0.07 


V12(ov) 


1.84 


15.12 


0.08 


14.73 


0.09 


14.74 


0.07 


V13 


3.20 


14.75 


0.1 1 


14.47 


0.11 


14.40 


0.08 


V14 


2.97 


14.97 


0.05 


14.56 


0.06 


14.54 


0.06 



Table 3: Mean IR magnitudes for Cepheids in NGC 2031. 

4. Light Curves and PL Relations 

4.1. Light Curves and Mean Magnitudes 

The ad opte d periods for NGC 1 86 6 are taken from Welch et all 
i 19911) and lWelch & Stetson! dl993h . For a few Cepheids, peri- 



ods have been re-determined by Musel laet all (|2006) and show 
no sign of significant changes between the two epochs. The pe- 
riods of the NGC 2031 Cepheids have been taken from lMateol 
d!992h and from our own list. We folded the data using the 
adopted periods obtaining light curves in all the three bands 
and the color curves in J - K. Light and color curves are shown 
in Fig. for NGC 1866, and in Fig.|2b]for NGC 2031. No 
shift has been applied to the phases in order to have the zero of 
the phase as first point, as is usually done, since this has been 
found to be unimportant for the determination of mean magni- 
tudes. Note that variables V5 and V7 in NGC 2031 have only 
one measurement and do not appear in the light curve plots. 
As can be seen in the figures, the sampling is far from optimal 
and, in some cases, the scatter is relatively large. Nevertheless, 
the periodic shape is almost always clearly defined, with the 
exception of a few outliers that usually have larger magnitude 
errors. For these, we checked the images to verify the pres- 
ence of glitches or problems and we found that always, in these 
cases, there was either a problem of crowding, especially in the 
J filter, combined with a larger FWHM and/or patches of cos- 
mic rays or dead pixels that affected the measurement. In fact, 
the points that are more scattered around the fiducial light curve 
are invariably from runs #1 and #3, when the Las Campanas lm 
IR camera had poor cosmetic properties (namely a large patch 
of bad pixels in the right part of the chip). The points with the 
largest errors have been removed while we kept the others for 
the determination of the mean megnitude. 



As expected, the global amplitude of the oscillations is of 
the order of a few tenths of magnitude. We calculated mean 
magnitudes by averaging the data, rather than fitting a template 
to the light curves. This method is preferred for its higher speed 
with respect to the classical, more accurate, template curve fit- 
ting, because the amplitudes of the Cepheids' light curves are 
relatively small and our data provide a sparse but relatively uni- 
form phase coverage. The different uncertainties from variable 
to variable are mostly due to the scatter along the light curve 
and to individual nonuniformities in period sampling. Objects 
having light curves with the poorest sampling suffer from a 
larger uncertainty in the mean magnitude. Moreover, since the 
light curve sampling is sparse, fitting to a template light curve 
does not reduce significantly the error on the mean magnitude. 
In order to reduce the effect of outliers and obtain a robust es- 
timate, an algorithm based on a biweight location estimate has 
been used. The biweight location is a robust estimator of the 
"mean" which uses weights defined as a function of the devia- 
tion from the median and calculated in an iterative way. In this 
context, robustness means that its sensitivity to single outliers 
is strongly reduced with respect to classic estim ators. For a de- 
tailed description of the method (see, e.g. Mo steller & Tukevl 
119771) . The average error on the mean magnitude is of the or- 
der of the error on the single measurement and slightly larger 
than that expected from a template fitting of a well sampled 
light curve. Table|2]records mean magnitudes in the three bands 
for 15 Cepheids of NGC 1866. Data for the 13 Cepheids in 
NGC 2031 are given in Tabled 
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4.2. Period-Luminosity Relation 

In order to construct PL relations in all the near-IR bands, 
we first need to adopt values for the cluster reddenings. The 
case of NGC 1866 has been widely discussed in the past. 
A value of E(B - V) _ 0.12 mag ha s recently been pro- 
posed in iGroeneweeen & Salarisl d2003l) from the We s enheit 
function, but this value was questioned by IStorm et alJ J2005) 
based on an IR surface brightness analysis. We therefore use 
the cano nical value of E( B - V) = 0.06 mag for this clus- 
ter (see IStorm et"alll2005l and references therein). It should 
be pointed out that, by adopting Aj = 0.9546 x E(B - V), 
A H = 0.6110 x E (B - V), A Ks = 0.3842 x E(B - V) (see 
ICardelh et al.ll 19891) . the difference in using one value or the 
other values is of the order of 0.04 mag in J, and 0.02 mag in 
Ks. By contrast, NGC 2031 lies in a crowded region close to 
the LMC bar, so its reddeni ng is corres pondingly higher. Here 
we adopt a value of lMould et aDdl993l) E(B - V) = 0.18 mag, 
leading to correction terms of: Aj = 0.17 mag, A H = 0.1 1 mag, 
A Ks = 0.07 mag. 

Given the small range of periods covered by our program 
Cepheids, a fit of the PL relation to these objects only is of 
limited use, due to the large uncertainties on the derived coeffi- 
cients. As an example, the PL fit in the K band for NGC 1886 
gives: M K (1866) = -2.79 + 0.31 * log(P) + 15.78 + 0.15. 

Our sample can instead be used for a compar i son w ith 
field Cepheids, namely the sample of IPersson et alJ d2004l) in 
which Cepheids with period shorter than ~3 days are absent. 
This comparison has two goals: (1) to check if the cluster 
and field Cepheids obey the same relation and, if so, to de- 
rive an overall PL relation extending over an enhanced period 
range; and (2) to check if the break in the slope of the PL 
relation at about P = lOd , discussed by Ngeo w et al.l J2005); 
Kan bur & Ngeowl d2006l) is present in our sample. Figure |4] 
shows PL relations for the sample discussed in this work, the 
IStorm et alJ d2005l) variables not in cluded in our sample, and 
the sample of Pers son et alJ d2004l) . Cluster Cepheids appear 
in the plot, as expe cted, as a extension of the sequence of 
IPersson et alJ d20Q4t) and seem to describe the same PL rela- 
tion, as was previously noted by Storm et al. (2005) (see their 
Fig. 15). 

Figure [5] is a magnification of the previous figure at 
log(P) < 1 which shows this fact more clearly. Note that two of 
the variables in the NGC 203 1 sample have longer period than 
the bulk of the cluster Cepheids and are likely field variables. 
If confirmed, this would provide further evidence that field and 
cluster Cephei ds describe the same rel ation. 

In addition, iKellerfe Wood! (2006) report a mean metallic - 
ity of [Fe/H] = -0.34 + 0.03 for their field sample, which is 
quite close to the metallicities of our two target clusters. Of 
course, the effect of metallicity on the PL relation — if such 
a dependence does indeed exist — is greatly reduced in the 
near-IR, especially in the K band. Thus, the fact that our clus- 
ter Cepheids extend the tight relation of IPersson et H]<|2004) 
toward short periods is perhaps not surprising. 

A new fit of the PL relation using the entire sample has 
been performe d in each of th e three IR bands, and the results 
compared with Per sson et alJ d2004l) . The output is shown in 



Fig.|5]as the red dashed line. The derived coeffic ients (see Fig. 
are essentially identical to the ones found by IPersson et alJ 
<2003>. 

There is an ongoing debate over possible non-linearitie s 
in the PL relation of L MC Cepheids. iNgeovv et al.1 d2005l) : 
iKanbur & Ngeowl d2006l) and ISandage et al.l d2004 reported 
that the PL relation of the LMC Cepheids shows evidence for 
a change in slope at about P = 10 days. This effect is most 
appar ent in the optical (i . e, V-, R- and I-band) PL relations, 
whilelN geow et al.ld2005h : lKanbur & Ngeowl (Eo06) show that, 
at least in th e J and H bands, it is present also in the IR. On the 
other hand. iGieren et al] |2005), by using IR surface bright- 
ness technique, found a constant slope, in agreement with the 
trend for Galactic Cepheids. From Figures |4] and |5] it can be 
seen that the slope looks essentially constant in all the three 
IR bands. In fact, by adding the short-period cluster Cepheids 
to the Persson sample, the slope is virtually unchanged, if not 
slightly flatter than the original Persson fit. Hence, the new fit 
with the extended period range seems to confirm the hypothesis 
of constant slope in the LMC Cepheids. 
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We have also made a comparison with the predicted PL re- 
lations from nonlinear convective pulsation models, for metal- 
licities and period ranges appropriate for o ur sample (see 
iBrocato et ai1l2004t iBono et aljfl999bl 120021 and references 
therein). The theoretical PL relations in J and K have similar 
values for the slope, 

Mj = -3.00 ±0.04 log P- 2.40 + 0.06 



smaller than the values predicted on the basis of the canonical 
evolutionary scenario; see, e.g., Brocato et al. 2004,Caputo et 
al. 2005), then the theoretical fundamental PL relations become 



M, 



-3.11 ±0.051ogP- 2.05 ±0.07 



M K = -3.21 + 0.04 log P - 2.37 ± 0.05, 



M K 



-3.20 ± 0.05 logP - 2.66 ± 0.05, 



with intrinsic dispersions of 0.15 mag and 0.13 mag, respec- 
tively. Comparison of theoretical and observed relations yields 
an estimate for the DM of LMC: (m - M) = 18.63 ± 0.11 mag 
for J and (m - M) = 18.62 ± 0.10 mag for K, that is slightly 
higher than the most recent estimates for the DM of the LMC 
obta ined from comparison wi t h theoretical models (1 8.54 mag 
iMarconi & Clementinil Eoolt iKeller & Wood! 120061) or wit h 
observation al methods (18.52 mag [Clementini et al. 2003), 
(18.56 mag lGieren et al.l2005l) . In particular. lAlvesI d2004l) has 
reviewed a series of distance measurements, finding an average 
DM of 18.50 ±0.02 mag. 
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Fig. 4: Overall P L relations for the field Cepheids of 
IPersson et alJ d2004l) and the two targ et clusters. The black dot- 
ted line represents the PL relation of IPersson et alJ J2004I) . the 
red dashed line shows the overall fit obt ained with the com - 
bined cluster sample and the sample of IPersson et alJ (|2004). 
The symbols for Cepheids are given in the plot. 



with an instrinsic dispersion of 0.19 mag and 0.15 mag. 
By applying these relations to the fundamental Cepheid 
data for NGC 1866 we obtain the mean apparent dis- 
tance moduli (m - M)o = 18.45 ± 0.1 1 mag in J and 
(m - M)o = 18.42 ± 0.10 mag for K, in better agreement 
with the short distance scale. Note that these are random 
errors that include all the observational and theoretical source 
of random uncertainty but do not take into account possible 
systematics due to physics used in the models. These include 
uncertainties arising from the treatment of convection (i.e., 
the mixing length parameter), the mass-luminosity relation 
and the atmospheric models used to obtain observed values 
from theoretical luminosity and temperature. The effect of 
the equation of state is negligible, and varying the opacity 
does not produce effects on the models, provide d that one 
uses the most recent compilations ((see IPetroni et al] 12003. 
for a discussion of this dependence). The effects of mixing 
length is under investigation (Marconi et al., in preparation) 
and preliminary results indicate that by increasing a from 1 .5 
to 1.8, the PL relation slope becomes slightly steeper. Such a 
change is expected to have a relatively minor impact on the PL 
relations derived above. 
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The zeropoints of theoretical PL relations depend on both 
the physical and numerical assumptions required for the con- 
struction of the pulsation models. If over the same period range, 
mildly overlurninou s pulsation models (see lBono et alj fl999b: 
IBrocato etai1 f2004) are adopted instead of the canonical ones 
(following recent suggestions that Cepheid masses should be 



logP 



Fig. 5: Magnified view of the PL relation for log(P) < 1. Lines 
and symbols are the same as in Fig.0] 
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5. Discussion and Summary 

We have presented IR light and color curves for 15 and 
13 Cepheids in the young LMC clusters NGC 1866 and 
NGC 2031, respectively. The light curves are of sufficient qual- 
ity to derive mean J, H, and Ks magnitudes, allowing us to com- 
pare our results on cluster Cephe ids with the an alogous study 
done on field Cepheids by Perss on et al.l (120041). with the oret- 
ical models bv lBono et ail dl999bl) : lBrocato et aD d2004 and, 
for NGC 1866, with previous studies from the literature. 

The good agreement we find wit h field Cepheids has al- 
lowed us to extend the PL relation of IPersson et alJ J2004I) to 
shorter periods. A redetermination of the PL relation yields 
slopes that are in agreement with recent theoretical models. 
To obtain agreement with recent determinations of the LMC 
DM that gi ve the same value (i.e., (m - M )p = 18.54) for both 
RR Lyrae jMarconi & C lementinil l2005h and with Cepheids 
jKelle r & Wood 2006[> would require overluminous models. In 
this case, whether the overluminosity is due to convective core 
overshooting or mass-loss is still a point of active debate. 

A comparison of results obtained on NGC 1866 by dif- 
ferent authors and different methods shows that some authors 
favour canonical models with the presence of binaries, while 
others conclude a cer tain amount of core overshooting is re- 
quired. For example, Walk er et alJ 11200 ll) use MS fitting on 
HST data to derive a DM of (m - M) = 18.35 + 0.06 by us- 
ing the observed Hyades sequence and theoretical models to 
take i nto account metallicit y differences. Groun d-based stud- 
ies bv lTesta et alJ i 19991) and Barmi naet alJ J2002I) used essen- 
tially the same data but yielded the opposite conclusion about 
the overshooting issue. The only point of accordance seem to 
be the presence of bina ries at the tip of the MS. A later rede- 
termination by Salar is et al. (2003) seems to confirm the DM 
of NGC 1866 found bv lWalker et al.1 fcOOlh and finds that the 
surrounding field has a longer DM of (m - M)o = 18.53 + 0.07 
(i.e., presumably due to the tilt of the LMC intermediate-age 
disk). 

In this work, the PL relations for the two clusters give es- 
sentially the same results, indicating that the two clusters have 
nearly the same DM. The structure of the young LMC — 
w here by "young" we m ean the sample of Cepheids analyzed 
in IPersson et alJ J2004) — may be descr ibed as a tilted disk, 
(see, e.g., the discussion in Alves 2004) with the north-east 
quad rant closer to t he Earth, and with some amount of thick- 
ness (Ivan der Marel et afll2002l) . The fact that the two clusters 
share essentially the same DM is consistent with NGC 1866 
lying behind the disk and at the same distance of the bar where 
NGC 203 1 is located. Although our data cannot solve this im- 
portant dilemma, it is worth emphasizing some outstanding is- 
sues that need to be resolved. Probably the most crucial point 
is the reddening — a better determination of this important pa- 
rameter is clearly needed. On the other hand, metallicity has a 
less dramatic effect, particularly at near-IR wavelengths. 

New models have recently made some progress towards 
resolving the longstanding discrepancy between evolutionary 
and pulsational masses; in fact, the same models are applied 
to cluster and field Cepheids. Finally, since the largest uncer- 
tainties seem to arise from the current observational material, 



the ultimate solution to the LMC distance problem will almost 
certainly require increased samples that include Cepheids in 
additional LMC clusters. 
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